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A  study  has  been  made  of  methods  of  solution  of  Maxwell's 


equations  in  inhomogeneous  Isotropic  plasmas.  Results  of  that  study  are 
given  here.  The  results  include  a  discussion  of  a  general  method  of 
reduction,  for  a  variety  of  coordinate  systems,  of  Maxwell's  equations  to 
scalar  differential  equations  for  the  separate  field  components.  In 
cartesian  coordinates,  reflection  and  transmission  coefficients  are  obtained 
for  plane  waves  propagating  obliquely  through  a  stratified  slab  of  plasma 
of  finite  thickness  and  through  a  semi-infinite  stratified  plasma.  Some 
discussion  is  included  concerning  total  reflection  and  Brewster's  Law  in 
connection  with  stratified  media. 
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MATHEMATICAL  METHODS  IN  THE  STUDY  OF 
WAVE  PROPAGATION  IN  INHOMOGENEOUS  MEDIA 


1. _ INTRODUCTION. 

The  material  to  be  discussed  in  this  report  has  been  developed  to 
deal  with  some  mathematical  aspects  of  the  interaction  of  electromagnetic 
waves  with  ionized  media.  It  is  concerned  primarily  with  methods  of  analysis 
of  Maxwell's  Equations  in  a  material  medium,  together  with  such  additional 
equations  as  may  be  needed  to  characterize  the  interactions  of  the  waves 
with  the  media,  with  the  aim  of  obtaining  from  the  set  of  equations  certain 
formulae  describing  the  way  in  which  the  waves  are  reflected  from,  absorbed 
in,  or  transmitted  through  the  medium.  In  addition,  questions  of  energy 
propagation  and  wave  front  analysis  are  considered.  The  emphasis  throughout 
is  on  mathematical  methods  for  describing  the  behavior  of  waves  in  media 
whose  electromagnetic  properties  vary  in  space.  The  methods,  though 
originally  developed  for  the  analysis  of  electromagnetic  wave  propagation 
in  inhomogeneous  plasmas,  may  equally  well  be  applied  to  inhomogeneous 
dielectrics  or  to  the  study  of  acoustical  wave  propagation  in  inhomogeneous 
fluids. 

The  system  of  equations  which  will  be  examined  in  this  report  will 
be  assumed  linear.  The  latter  assumption  is  a  drastic  one.  However, 
inasmuch  as  the  solution  of  systems  of  nonlinear  equations,  even  in  the 
case  of  homogeneous  media,  involves  the  use  of  analytical  tools  poorly 
understood  even  by  specialists  in  the  field  of  nonlinear  phenomena,  it 
seems  wise  to  tackle  first  the  problems  associated  with  linear  equations 
in  inhomogeneous  media.  The  latter  sometimes  present  knotty  problems 
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with  which  we  hope  to  deal  in  this  report;  but  the  analytical  tools  (e.g., 

WKB  methods)  available  for  their  treatment  are  much  better  understood, 

The  introduction  of  Maxwell's  Equations  into  the  analysis;  and 

their  subsequent  manipulation;  leads  to  a  vector  differential  equation  for 

★ 

the  electric  field  E.  The  latter  equation  is  equivalent  to  a  Bet  of  threa 
scalar  equations;  with  variable  coefficients;  in  which  the  components  of 
E  are  "coupled."  The  "uncoupling"  of  these  field  components,  that  is, 
the  derivation  of  scalar  differential  equations  for  the  individual  field 
components,  is,  in  general,  difficult,  especially  when  one  is  dealing 
with  coordinate  systems  other  than  cartesian.  In  Appendix  III  we  have 
given  a  procedure  for  "uncoupling"  the  field  components.  The  procedure 
is  certainly  not  applicable  to  the  general  case  but  assumed  that  the  medium 
in  question  varies  as  a  function  of  only  one  space  coordinate.  Also,  the 
orthogonal  curvilinear  coordinate  system  to  which  the  procedure  applies  is 
not  the  most  general.  However,  in  spite  of  these  limitations,  the  procedure 
seems  to  be  of  sufficiently  wide  application  to  merit  consideration. 

In  the  treatment  of  reflection  and  transmission  coefficients  a 
difficulty  arises  in  the  case  of  spatially  non-constant  media  in  that  the 
prescribing  of  the  solutions  of  the  relevant  differential  equations  may 
require,  when  one  is  discussing  a  semi-infinite  medium,  the  prescribing 
of  boundary  conditions  at  infinity,  conditions  which  are  not  always  easily 
justifiable.  A  method  of  approach  to  this  difficulty  of  prescribing 
boundary  conditions  at  infinity  is  presented  in  Appendix  V.  While  this 

i 

*  Capital  letters  without  subscripts  will  be  used  to  represent  vectors  in 
this  report.  Otherwise,  all  symbols  refer  to  scalars. 
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method  also  presents  difficulties  of  interpretation,  it  appears  to  shed 
some  light  on  the  source  of  the  boundary-condition  difficulty  Just 
mentioned  and  it  also  appears  to  be  superior  to  any  method  presented  in 
the  literature  on  wave-propagation  analysis  in  inhomogeneous  media. 

II. _ DISCUSSION  OF  THE  FUNDAMENTAL  EQUATIONS. 

In  the  subsequent  discussion,  MKS  units  will  be  used  throughout. 

The  behavior  of  the  electromagnetic  disturbance  in  the  medium  under 

consideration  will  be  described  mathematically  by  Maxwell's  Equations.  The 

ionized  medium  under  consideration  will  be  considered  to  be  a  dilute  plasma 

at  low  pressures  and  isotropic,  so  that,  among  other  things,  polarizability 

effects  can  be  neglected,  and  the  dielectric  permittivity  and  magnetic 

permeability  of  the  medium  will  have,  to  a  good  approxmation,  the  values 

which  those  quantities  have  in  the  vacuum,  and  these  values  will  be  denoted 

by  €q  and  respectively.  Then  it  will  be  possible  to  write  down  the 

so-called  "constitutive  relations"  D  -  e  E  and  B  ■  u  H.  Hence,  the  only 

o  o 

way  in  which  the  influence  of  the  material  medium  is  introduced  into  the 
problem  is  through  the  current  density  J  appearing  in  one  of  Maxwell's 
Equations.  The  introduction  of  the  constitutive  relations  into  Maxwell's 
Equations  leads  to  the  pair  of  equations 

V  *  H  -  «0  "ft  +  J  (a) 

U) 

V  x  E  -  -  nQ  -§f  .  (b) 

The  field  variable  H  can  be  eliminated  from  the  latter  pair  of  equations 
by  the  familiar  procedure  of  taking  the  curl  of  the  second  equation  and 
replacing  the  resulting  quantity  Vx  H  by  its  equivalent  from  the  first 
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equation.  The  equation  resulting  from  the  elimination  of  H  is 


0)(o>  +  iv) 


and  Equation  (6)  then  becomes 

Vx  Vx  E  ■  kE-  (8) 

The  coefficient  k  in  Equation  (8)  is  now,  in  general  to  be 
considered  a  complex  function  of  the  space  coordinates.  When  k  is  constant, 
or  nearly  constant  the  solution  to  (8)  can  be  discussed  in  terms  of  the 
"geometrical  optical"  approximation  in  which  wave  fronts  and  rays  can  be 
described  through  the  use  of  Hamilton's  characteristic  function  S  (eikonal 
of  Bruns),  discussed  in  Appendix  II.  The  quantity  k  plays  the  role  of  a 
spatially  varying  index  of  refraction.  However,  as  k  gets  close  to  zero  the 
geometrical  optical  approximation  becomes  less  and  less  accurate,  and  some 
new  method  of  approximation  must  be  used.  In  general,  however,  the  problem 
of  ascertaining  the  behavior  of  an  electromagnetic  wave  in  an  inhomogeneous 
medium  requires  a  detailed  solution  of  the  vector  wave  equation  (8),  the 
features  of  the  solution  depending  strongly  on  the  prescribed  functional 
form  of  k.  A  more  detailed  discussion  of  the  implications  of  various 
approximate  solutions  can  be  found  in  reference  2. 


III.  PRELIMINARY  DISCUSSION  OF  THE  PROCEDURE  IN  SOLVING  THE 

FUNDAMENTAL  EQUATION. 


As  yet,  no  commitment  has  been  made  in  the  above  discussion  to  the 
following,  which  will  be  discussed  in  the  order  given: 

A.  The  restrictions  needed  on  the  functional  form  of  k. 

B.  The  types  of  waves  (plane,  cylindrical,  spherical,  etc.) 
in  terms  of  which  the  solutions  to  Equation  (8)  may  be 
constructed. 
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C.  The  shapes  of  the  boundaries  of  the  media  with  which 

c 

the  waves  interact.  Media  may  be  assumed  infinite  in 
one  or  more  directions. 

D.  The  choice  of  coordinate  systems  in  terms  of  which  the 
solutions  to  (8)  can  be  expressed. 

E.  The  boundary  conditions  which  must  be  satisfied  by  the 
solutions  to  (8). 

Each  of  the  preceding  five  topics  will  now  be  discussed  in  turn. 

A. _ The  Functional  Form  of  k. 

Equation  (8)  is  a  vector  equation  for  the  unknown  electric 
field  distribution  E.  As  such,  the  equation  is  equivalent  to  a  set  of 
three  simultaneous  linear  scalar  partial  differential  equations  for 
three  unknown  variables,  the  mutually  orthogonal  components  of  E,  relative 
to  a  prescribed  three-dimensional  orthogonal  coordinate  system.  In  general 
the  components  of  E  are  said  to  be  "coupled"  in  each  of  the  three  afore¬ 
mentioned  scalar  equations.  In  Appendix  III,  a  discussion  is  given  of  the 
system  in  special  cases. 

The  following  observations  should  be  made,  however,  concerning 
the  system  III(c)  presented  in  Appendix  III. 

1.  Solutions  of  III(c)  which  are  amenable  to  further 
manipulation  for  the  mathematical  derivation  of  such  "observables"  as 
phase  shift,  absorption,  reflection  coefficients,  transmission  coefficients 
etc.,  generally  require  that  two  of  the  three  unknown  field  components  be 
eliminated  from  the  system  III(c),  yielding  a  single  scalar  linear  homo¬ 
geneous  partial  differential  equation  for  a  single  field  component. 


Alternatively,  one  may  introduce  new  field  variables,  usually  called  "po¬ 
tentials"  (Hertz  potential,  scalar  and  vector  potentials)  and  obtain  scalar 
equations  for  each  of  them. 

2.  The  carrying  out  of  the  aforementioned  elimination  is 
impossible,  in  general,  when  the  coordinate  system  and  the  functional 
form  of  k  are  arbitrarily  prescribed.  In  fact,  the  elimination  appears, 
in  the  present  state  of  the  art,  to  be  feasible,  in  general,  only  in  the 
case  in  which  k  is  a  function  of  only  one  coordinate. 

3.  If  the  dependence  of  k  and  of  E  on  the  coordinates  is  not 
suitably  chosen  to  conform  to  *-he  shapes  of  the  bounding  surfaces  of  the 
media  in  question,  a  further  difficulty  arises  in  setting  up  and  solving 
the  relations  engendered  by. the  boundary  conditions  at  the  bounding 
surfaces. 

For  the  reasons  just  given,  the  only  inhomogeneous  medium 
problems  which  have  so  far  yielded  substantially  to  mathematical  and 
numerical  treatments  are  those  involving  so-called  "stratified  media." 

The  latter  media  may  be  regarded  as  being  built  up  of  "lamina"  of  homogeneous 
material,  or  they  may  be  continuously  varying  in  one  coordinate  direction, 

4.  In  discussions  of  homogeneous  media  it  is  frequently  desirable 
to  resolve  the  solutions  of  (8)  into  "outgoing"  and  "incoming"  waves.  In 
the  very  simple  situation  involving  two  homogeneous  media  separated  by  an 
infinite  plane  interface,  the  incoming  wave  in  one  of  the  media  is  often 
discarded  completely,  the  physical  reason  being  that  no  source  of  energy 
exists  at  infinity  in  that  medium.  However,  even  in  the  simplest  case  of 
waves  propagating  in  the  direction  of  variation  of  a  medium  varying  in 
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only  one  direction,  it  is  not  always  possible  to  express  the  solution  of 
(8)  unambiguously  in  terms  of  two  solutions,  one  of  which  can  be  identified 
as  an  "outgoing"  disturbance,  while  the  other  can  be  regarded  as  an  "incom¬ 
ing"  disturbance.  Some  examples  illustrating  this  point  are  cited  in  Appendix 
IV. 

J3. _ Types  of  Waves  Used  in  Setting  Up  Solutions. 

The  first  and  most  important  remark  to  be  made  here  is  that 
a  wave  of  one  type  can  be  represented  as  a  (usually  infinite)  sum  of 
waves  of  another  type.  Specifically,  a  spherical  wave  can  be  represented 
as  the  sum  of  plane  waves  (cf.  ref.  3,  p.  238  and  ref.  6,  chap.  VI).  Ihis 
amounts  to  writing  the  expression  for  the  spherical  wave  as  a  Fourier 
series  or  as  a  Fourier  integral.  Similar  remarks  can  be  made  for  other 
types  of  waves.  The  selection  of  the  type  of  wave  in  terms  of  which  the 
solution  to  the  wave  equation  will  be  expressed  will  generally  be  determined 
by  shapes  of  the  boundaries  at  which  boundary  conditions  are  prescribed. 

Thus  cylindrical  boundaries  will,  in  general,  lead  to  the  choice  of  cylindri¬ 
cal  waves  as  the  most  convenient  "basis"  for  the  solution. 

£. _ Boundary  Shapes. 

The  vast  majority  of  boundaries  encountered  in  wave  problems 
are  either  plane,  cylindrical,  or  spherical.  Many  problems  of  practical 
interest  can  be  framed,  perhaps  as  an  approximation,  in  such  a  way  that 
the  boundaries  in  question  are  plane.  Furthermore,  while  problems 
involving  other  boundary  shapes  can  be  solved,  their  solution  involves 
the  use  of  more  sophisticated  analysis  (notably  intricate  contour  inte¬ 
grations  in  the  complex  plane)  which  would  tend  to  bury  the  physical 
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conclusions  in  a  welter  of  analytical  detail.  Therefore,  the  surfaces  at 
which  boundary  conditions  will  be  prescribed  in  this  report  will  always  be 
taken  to  be  plane. 

I). _ Choice  of  Coordinate  System, 

It  is  a  familiar  fact  that  the  solution  to  boundary  value  problems 
is  often  facilitated  when  the  coordinate  system  is  chosen  so  that  the 
boundary  surfaces  are  coordinate  surfaces.  Since  this  report  is  committed 
only  to  the  discussion  of  boundary  value  problems  involving  plane  boundaries, 
all  boundary  value  problems  in  the  report  will  be  treated  in  terms  of 
cartesian  coordinates.  However,  some  attention  will  be  given,  especially 
in  Appendix  III,  to  the  manipulation  of  Equation  (8)  in  other  coordinate 
systems. 

IS. _ Boundary  Conditions. 

Following  the  conclusions  arrived  at  in  Stratton  (ref.  6,  p.  37), 
we  shall  assume  that  the  tangential  components  of  E  and  H  are  continuous 
across  any  boundary  at  which  the  conductivity  is  finite.  The  continuity 
conditions  will  be  used  in  Appendix  V,  where  reflection  and  transmission 
coefficients  are  derived  for  plane  waves.  More  difficult  decisions  are 
posed  in  considering  the  behavior  of  the  solutions  to  the  wave  equation 
at  infinity  in  cases  in  which  the  inhomogeneous  medium  extends  to  infinity 
in  the  direction  of  the  gradient  of  k.  In  view  of  the  conclusions  implied 
by  the  discussion  given  in  Appendix  IV,  only  those  cases  are  considered  in 
which  the  solutions  of  the  wave  equation  have,  for  large  values  of  the 
argument,  asymptotic  forms  which  have  a  wave-like  character. 
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IV. 


FURTHER  DISCUSSION  OF  THE  SOLUTION  OF  THE  FUNDAMENTAL  EQUATION. 
REDUCTION  TO  EQUATIONS  FOR  EACH  FIELD  COMPONENT. _ 


In  Appendix  III  a  scalar  partial  differential  equation  has  been 
derived,  under  certain  restrictive  assumptions,  for  each  of  the  components 
of  the  electric  field  E.  The  results  are  summarized  below.  For  each 
coordinate  system  the  scale  factors  (described  in  Appendix  III)  are  given, 
and  the  assumptions  are  given  under  which  the  differential  equations  for 
the  various  field  components  were  derived.  The  Differential  equations  are 
special  cases  of  III(j),  III(k),  and  III(p). 

A.  Cartesian  Coordinates  (x,y,z).  Scale  factors:  h  *  h  *  h  *  1. 

_  ' _ x  y  z 

The  quantity  k  is  assumed  to  be  a  function  of  x  only.  The  field 

components  E  ,  E  ,  E  are  assumed  to  be  independent  of  y.  They  satisfy 
x  y  z 

the  following  equations: 


(D  2  +  D  2  +  k)  E  =0 
x  z  y 

D  (  f  D  )  (  kE  )  +  D2E  +kE  *0 

X  K  X  X  Z  X  X 


<Dg  +  k)  )  k 


2  2 
D  +  D  +  k 
x  z 


} 


+  ’<Dz  W  (Dxk> 


(a) 

(b)  (9) 


(c) 


B.  Circular  Cylindrical  Coordinates  (r,  9,  z).  Scale  Factors: 

_ hr  "  hz  '  1>  h8  *  r- _ 

The  quantity  k  is  assumed  to  be  a  function  of  r  only.  The  field 


components  E  ,  E  ,  E  are  assumed  to  be  independent  of  8.  They  satisfy  the 
r  o  z 

following  equations: 


10  - 


1 


(D  -D+-D2+-)  (rEa)  -  0 
'  r  r  r  r  z  r  0 


(a) 


D  (  r  D  )  (rkE  )  -  -  b  D  (rkE  )  +  r  (D  2  +  k)E  -  0 

r  k  r'  '  r'  r  |_kr  r'J  z  '  r 

[k  [V°r  +  rDz2  +  k]  ) 


(Dz  +fc) 


E  + 
z 


(b) 

(10) 


+  r(D  D  E  )  (D  k)  -  0 
v  z  r  z'  '  r  ' 


(c) 


C. 


Spherical  Coordinates  (r,  Cp  ,  0).  Scale  Factors:  h  -l, 
_ -  r  aln  8,  hfl  »  r.  _ * _ 

The  quantity  k  is  assumed  ;to  be  a  function  of  r  only.  The  field 
components  E^,  E  ^  ,  Eg  are  assumed  to  be  Independent  of  4>  .  They  satisfy 
the  following  equations: 


sin  0  D  +  D. 

r  0 


r  sin  .0 


Dg  |  +  k  sin  0  ^  (r  sin  OE^  )  -  0  (a) 


sin  9Dr(^Dr)  (r2kEr) 


0  sin  0  0 


D.  +  r  k  sin  9 


0  (b) 

(11) 


D„  sin  9  D  E.  =  (Da  sin  0DQ  +  k  sin  0)  E 
y  r  y  y  w  r 


(c) 


The  latter  equation  is  from  Equation  Ill(i^)  of  Appendix  III. 

The  sets  of  Equations  (9),  (10),  and  (11)  are  suitable  for  appli¬ 
cation  in  a  variety  of  interior  and  exterior  boundary  value  problems,  e.g., 
waveguide  and  scattering  problems.  The  limitation  that  the  field  components 
must  be  independent  of  u^,  however,  restricts  their  use  somewhat.  As  was 
mentioned  in  Appendix  III,  analogous  sets  of  equations  can  be  derived  with¬ 
out  the  aforementioned  restriction,  but  the  equations  can  become  quite 
cumbersome.  Only  the  set  (9)  will  receive  further  attention  in  this  report. 
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V. 


THE  REFLECTION  AND  TRANSMISSION  COEFFICIENTS  FOR  PLANE  WAVES. 


In  this  section  formulas  will  be  given  for  the  reflection  and 
transmission  coefficients  for  a  plane  wave  obliquely  incident  on  the  plane 
interface  separating  a  homogeneous  medium  from  a  stratified  medium.  The 
formulas  will  be  derived  in  Appendix  V  on  the  basis  of  the  following 
assumptions: 

(i)  The  medium  is  stratified  in  the  x-direction;  i.e., 
k  is  a  function  of  x  alone. 

(ii)  The  stratified  medium  is  bounded  by  infinite  parallel 
planes  at  x  ■  0  and  at  x  *  d.  The  plane  wave  is 
assumed  to  be  incident  on  the  x  *  0  plane,  from  the 
negative  x-direction. 

(iii)  The  plane  of  incidence  of  the  incident  plane  wave  is 
the  x-z  plane  so  that  all  field  components  are 
independent  of  the  y-coordinate  and  can  therefore 
be  determined  from  the  system  (9).  (There  is  no 
loss  of  generality  in  this  assumption,  since  the 
y-axls  can  always  be  taken  to  lie  in  the  x  »  0 
plane  and  orthogonal  to  the  direction  of  propagation 
of  the  plane  wave.) 

(iv)  The  electric  and  magnetic  field  components  E^,  Ez,  H^, 

Hz  are  continuous  everywhere.  The  boundary  conditions 
on  these  field  variables  are  that  they  are  continuous  at 
x  -  0,  and  at  x  ■  d,  for  every  value  of  z. 

The  derivation  of  the  reflection  and  transmission  coefficients 

will  employ  the  traditional  method  of  considering  two  separate  cases: 

the  case  in  which  the  incident  electric  field  is  polarized  perpendicular 

* 

to  the  plane  of  incidence  (TE-wave  case)  and  the  case  in  which  the  magnetic 
*  These  designations  TE-wave  and  TM-wave  are  adopted  from  Born  and  Wolf. 
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field  is  polarized  perpendicular  to  the  plane  of  incidence  (IM-wave  case). 
Since  the  incident  electromagnetic  field  can  always  be  resolved  into 
components  described  by  these  two  cases,  the  reflection  and  transmission 
coefficients  will  be  completely  determined  by  the  formulas  derived  in 
these  two  cases. 

A  sunmary  of  the  formulas  derived  in  Appe'ndix  V  is  given  below 
(refer  to  Figs.  1  and  2  in  Appendix  V  for  the  meanihg  of  the  symbols). 


=  ik^cos  0^ 


D2  =  ik3  cos  03 


TE- 

■WAVE  CASE 

4>+«» 

4)  (°> 

b\o) 

'  + 

f 

<J>  (o) 

A  <*> 
> 

(f>  (d) 

(fc'(d) 

+ 

<f>'  (d) 

p  (o) 

1  + 

p  (o) 

t)  <o) 

r  + 

<p  (o) 

Km 

( \ )  (d) 

1 

(P  (d) 

+ 

1 

P  (d) 

A  =  D3  +  (ik3  cos  0^  D2 
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-2ik^  cos  0^ 


•  i 


£+(d)  ^.(d)  -  4>  (d)  4>_(d)  -ik_d  cos  0 

~  _  A  -  ^ 


Al6 


*1 


P3  C08  ©3 

k^  cos  0^ 


G  ,  Mil 

**1  k3  cos  B3 


G  =  Mil  cog  e 

^2  k3  cos  H3 


TM-WAVE 

CASE 

1 

fi  (o) 

+ 

1 

/°  (o) 

/>  (d) 

+ 

/>  (d) 

'  /Mo) 

f(o) 

4(d) 

f>  (d) 

/°  (o) 
+ 


/*(d) 
'  + 


^  (o) 


(d) 


/7+(°)  /t(0) 


Z5  (d)  (d) 


Gl  + 


iSk^  cos  ei 


-Gi  +  <  “s  ei  > 


a 
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CL 


-2  cos  8 
“l  1 


t°.W  -  cos  9  /".(d) 


a. 


a. 


cos  0, 


(d)  -  -^141  cos  03  s*  (d) 


L  + 


+  J 


°L 


,2  cos  e 

*1 _ L 


/*(<*)  /°  (d)  -  ^  <d)  ^  (d) 


“lk3d  cos  0. 

OLe  3 


V11 


“l 
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a. 


at 


cos  0, 
cos  0, 


An  important  difference  between  the  two  cases  is  the  following: 

The  TE-wave  case  depends  on  the  form  of  the  functions  <^>+(x) 

and  (J)  (x)  which  are  solutions  to  the  differential  equation 

<^>  "  +  (k  -  y  )  <t*  =  while  the  TM-wave  case  depends  on  the 

form  of  the  functions /-’  (x)  and  (x)  which  are  solutions  to 

+  1  ”  2 

the  differential  equation  +  *  0  (cf. 

equations  V(a),  V(p),  V  (q)).  Note  that  the  two  differential 
equations  are  the  same  when  k  is  constant.  Hence,  this 
distinction  between  the  two  cases  is  introduced  trfien  we  allow 
k  to  vary  in  space. 

The  case  in  which  medium  II  is  homogeneous,  with  propagation  constant  k 2,  is 
so  instructive  for  its  lessons  in  the  inhomogeneous  case  that  we  introduce 
it  here.  For  this  case,  the  following  relations  hold: 

+  ik  x  cos  0„ 

(j>  (x)  *  e  2  =  ^  (x);  k^  sin  01  =  k2  sin  ©2  =  k3  sin  03* 
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TE-WAVE  CASE 


■  2k 2  cos  ©2  _k3  cos  ®3  cos  0^  C0S  +  ^2  008  ®2  sj['n  ^2d  COS  ®2^j 

D2  ■  21  ik^  cos  03  sin  (k^d  cos  ©2)  +  k2  cos  02  cos  (k^d  cos  @2> 

L  r  1  J 

2k ^  cos  0^  [^k3  cos  9^  -  k2  cos  0^ 


Now  note  that  for  example,,  is  the  amplitude  o£  the  wave  propagating 
backward  in  medium  II,  essentially  the  reflected  wave  from  the  interface 
between  medium  II  and  medium  III.  This  reflection  can  be  extinguished  by 
having  k^  =  k2,  whence  k^  cos  0^  =  k2  cos  9^  (consequence  of  Snell's  Law) 
and  B2  =  0.  The  latter  observation  suggests  a  means  for  our  determining 
the  behavior  of  the  disturbance  in  medium  II  when  the  latter  is  inhomo¬ 
geneous  and  semi-infinite  (d  — *  oo).  We  proceed  as  follows:  let  medium 


III  be  homogeneous  with  the  propagation  constant  r  k(d)  -  y  which  medium 
II  would  have  were  it  suddenly  to  become  homogeneous  at  x  *  d,  keeping  the 
value  of  the  dielectric  permittivity  which  it  attains  at  x  =  d,  as  dia¬ 
grammed  below 
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Then,  instead  of  characterizing  the  obliquely  propagating  wave  by  a  propa¬ 
gation  constant  cos  0^  in  medium  III,  we  replace  it  by  a  quantity  <7- (d), 
dependent  only  on  d.  It  can  be  seen  then  that  (reflection  amplitude  in 


medium  II) 


2ik1  cos  0X  [<£>+(d)  -  io-(d)  4>+(d)] 


Now,  suppose  that  the  ordinary  differential  equation  determining  <p  (x) 
admits  of  two  solutions  (f>  +(x)  and  <|>_(x)  with  the  property  that 


d^oo  (j)  (d)  =  +  <3-(d)  4>  (d). 

-  ±  J  ± 


That  is,  for  large  values  of  d,  we  may  replace  <j>  (d)  by  +  i  cr(d)  4>  (d), 

+  + 

to  a  good  approximation.  Then  %2~~*  ®  a8  d  — *  ®  •  Let  ub  now  calculate 
the  reflection  coefficient  in  medium  I,  for  large  d,  using  the  aforementioned 
assumption.  According  to  the  table  on  page  13 

1 

-D^  +  ik^  cos 

®1D2  '  4+(o)  +  ikl  cos  4>+( o) 

B1  D.  +  ik.  cos  0,D_  ^  Q  l  ,  N  Al* 

11  12  A  (o)  +  ik^  cos  0-  <p  (o) 

V+  1  1  + 


It  follows  from  the  latter  relation  that  *  |  ~  j  will  be  equal  to 

unity  (total  reflection  case)  if  (j)  (o)  =  0  or  if  (j)  (o)  *  0  or  if  cos 

+  .  + 

6.  =  0  (grazing  incidence).  When  medium  II  is  a  homogeneous  medium  so 
ik^x  cos  02 

that  (p  +(x)  *  e  ,  the  non-trivial  total  reflection  condition 

becomes  cos  @2  =  0,  in  agreement  with  the  familiar  result.  Note  that, 

from  the  continuity  conditions,  Snell's  Law  yields  k^  sin  0^  *  y  *^k(o)  -  ^-(o), 

The  condition  (o)  =  0  is  a  kind  of  eigenvalue  condition  which  will 
T  + 


17 


determine  ®-(o),  which  will  in  turn  determine  the  angle  0.  for  total 

i  i 

reflection.  Note,  also,  that  the  condition  4>  +(o)  ■  0  might  relate  the 
reflection  to  the  region  in  the  dielectric  where  the  index  of  refraction 
vanishes,  but  only  when  <$)  ,  (x)  is  approximately  representable  in  complex 
exponential  form  (propagating  wave)  near  x  =  0.  Otherwise,  almost  any 
condition,  unrelated  to  the  index  of  refraction,  might  arise.  All  of 
these  observations  apply  only  to  the  lossless  case,  of  course. 

The  above  remarks  have  all  applied  to  the  TE-wave  case.  In  the 


TM-wave  case,  it  is  possible  to  deduce  for  a  semi-infinite  medium  a 
condition  on  the  angle  of  incidence  0^  such  that  »  0.  The  specific 
angle  of  incidence  in  this  case  is  called  the  "Brewster  angle."  In  this 
case  we  must  argue,  not  on  the  functional  behavior  of  (j)  (x),  but  on  the 
behavior  of  a  function  /P  (x)  which  is,  in  general,  entirely  different  from 
cj)(x).  The  argument  here,  however,  parallels  the  preceding  one  and  yields 


the  condition 


/>  (o) 

*£*  CCS  8.  -  -1  -±j-r 

kl  1 

for  the  Brewster  angle.  When  /°+(x) 


ik2X  cos  0£ 

e  ,  the  latter  condition 


becomes  k2  cos  0^  =  k^  cos  @2*  Combining  this  with  Snell's  Law,  k^  sin 


@1  -  k2  sin  &2>  we  8et  sin  20^^  ■  sin  202*  In  the  non-trivial  case,  the 
latter  condition  implies  ©2  *  j  -  0^.  Substituting  this  in  Snell's  Law, 
we  get  k^  sin  0^  «  kj  cos  0^,  or  tan  0^  *  k2/k^,  yielding  Brewster's 


angle  for  the  case  in  which  medium  II  is  homogeneous  and  extends  to 


x  *  +  co  . 
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It  is  of  considerable  interest  to  note  that,  when  medium  II  is 


homogeneous  and  semi-infinite,  no  zero-reflection  situation  can  arise 
in  the  TE-wave  case  and  no  total  reflection  can  arise  in  the  TM-wave 
case.  However,  when  medium  II  is  allowed  to  be  inhomogeneous,  both 
reflection  situations  can  exist  in  either  case.  The  examples  are  too 
complicated  to  cite  here,  but  some  information  on  the  subject  can  be 
found  in  Reference  7,  page  19. 

VI. _ CLOSED-FORM  SOLUTIONS  OF  THE  WAVE  EQUATIONS. 

The  question  will  arise  as  to  what  functional  forms  for  k(x) 
will  permit  the  closed-form  solution  of  Equations  (9-a)  and  (9-b)  simul¬ 
taneously.  Aside  from  the  trivial  case  of  k(x)  constant  there  is  the  case 

Ox 

of  k(x)  proportional  to  e  where  C  is  any  constant.  The  solutions  in 
this  case  will  be  given  in  the  form  of  Bessel  functions.  In  general, 
exponential  variations  for  k(x)  seem  to  show  the  most  promise  in  this 
connection.  When  k(x)  is  a  linear  function  of  x.  Equation  (9-a)  admits 
closed-form  solutions  in  terms  of  Airy  or  Bessel  functions,  but  Equation 
(9-b)  does  not  seem,  in  this  instance,  to  admit  closed-form  solutions  in 
terms  of  known  functions.  The  whole  question  discussed  in  this  section  is 
still  under  investigation. 
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APPENDIX  I 
•* 

Dependence  of  J  on  E 


The  equation  for  the  current  density  which  will  be  used  in  this 
report  is  that  derived  by  Spitzer  (ref.  1,  p  21,  equation  2-12).  In  our 
notation  (J  *  current  density)  and  in  our  system  of  units  (MKS),  the  afore¬ 
mentioned  equation  of  Spitzer  is  as  follows: 


d  J  Z  ,  ■*  „  .  1 

-j-  =  E  +  v  x  B  + -  VP„  " 

dt  en  e 

e 


The  following  assumptions  will  be  made: 


x  b  -  a  J. 


i.  v  -  0  since  the  medium  is  assumed  to  be  at  rest  over  all. 

ii.  A  new  quantity,  <0^,  the  plasma  frequency,  defined  by 


-,  will  be  introduced. 


iii.  The  coefficient  ^  >  used  in  Equation  (la)  is  essentially 

a  measure  of  energy  dissipation.  We  introduce  in  its  place 

n  e^ 

the  dissipation  parameter  v  defined  by  v  ■  -  H  . 

m 

e 

iv.  The  most  drastic  assumption  is  that  the  term 
1  -* 

-  J  x  B  is  negligible  with  respect  to  the  other 

e 

quantities  appearing  in  Equation  (la).  This  assumption 
is  essentially  that  the  field  amplitudes  are  small 
(small-signed  theory)  and  that  the  medium  is  isotropic. 


Equation  (la)  then  reduces  to 


77  +  vJ  =  a)  2e  E  + 
dt  p  o 
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At  this  point  it  becomes  desirable,  mostly  for  the  purpose  of 

simplifying  the  mathematical  treatment,  to  assume  that  all  of  the  field 

quantities  vary  in  time  at  a  single  frequency  co,  so  that  for  any  field 

d  F 

variable  F  whatsoever,  ■  -icoF.  It  follows  from  the  latter  assumption 
that  Equation  (lb)  can  be  expressed  as  follows: 


(-ico  +  v)J  ■  co  e  E  + 
P  o 


(t) 


VP. 


or 


,  2 

r  i  i 

ico  e 

m 

P  o 

U  i 

e 

£1  + 

_  CO  +  iv  _ 

co  +  iv_ 

VP, 


(Ic) 


On  the  basis  of  elementary  kinetic  theory,  one  concludes  that  p  *  nkT 

©  6 

2 

or  pe  *  o)p  (eo^/e)T.  It  follows  that  <7  pe  may  make  an  important  contri¬ 
bution  to  J  whenever  ne  sufficiently  large.  The  inclusion  of  V  Pe 

in  J  would,  when  substituted  into  Equation  (4),  yield  an  inhomogeneous 
partial  differential  equation  for  E,  an  equation  with  quite  important 
consequences.  However,  in  any  event  it  would  be  necessary  to  solve  the 
homogeneous  equation  first  in  order,  for  example,  to  obtain  the  Green's 
Function  for  the  boundary-value  problem  in  question.  Hence,  it  has  been 
considered  desirable  here  to  consider  first  the  consequences  of  having 
y  Pe  *  0.  The  latter  assumption,  in  our  context,  is  tantamount  to 
considering  a  zero-temperature  plasma.  The  final  formula  for  J  which  will 

be  used  in  connection  with  Equation  (4)  is 

r  ,  2  ” 

iCO  € 

P.  o. 


co  +  iv 


E. 


*  cf.  e.g.,  A  Sommerfeld,  Thermodynamic s .  etc.  pp.  8-12. 
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APPENDIX  II 


Hamilton's  Characteristic  Function  (the  Eikonal  of  Bruns) 


A  brief  discussion  will  be  given  here  of  the  eikonal  (characterls- 

•ff 

tic  function  of  Hamilton)  defined  by  H.  Bruns.  Equation  (8),  Introduced 
in  the  main  body  of  the  report,  is 

Vx  ^7x  E  *  kE,  (Ila) 

where  k  is  a  scalar  function  of  the  coordinates  and  E  is  a  vector  function 
of  the  coordinates.  Taking  the  divergence  of  both  sides  of  (Ila)  yields 

0  -  V  •  (kE)  -  E  •  V*  +  k  V  •  E  (lib) 

whence 


V  *  E  -  -  E  •  ( V  k/k). 


(lie) 


Now,  with  the  help  of  (lie),  and  the  relation  Vx  Vx  E  *  ^7 (  V  *  E)  -  V  n, 
we  may  rewrite  (Ila)  as  follows: 

,2. 


V‘E  +  kE  =  -  V  E  •  — 1  . 


(lid) 


If  k  were  constant,  Equation  (lid)  would  reduce  to  the  familiar  vector 
Helmholtz  equation  for  E,  for  which  we  could  give  the  solution 
E  *  Ae^klx  +  k2^  +  k3z^  in  cartesian  coordinates,  where  k^2  +  +  k^2  =*  k, 

and  A  is  a  constant  amplitude.  With  the  latter  remark  in  mind,  we  try  to 
"fit"  a  solution  of  the  form 

E  =  AeiS  (lie) 

to  the  Equation  (lid),  where  A  and  S  are  now  permitted  to  vary  in  space. 

We  substitute  (lie)  into  (lid)  and  ascertain  the  conditions  under  which 
(lie)  will  be,  approximately,  a  solution  to  Equation  (lid).  The  solution 
(lie)  is  the  "geometrical  optical"  approximation  to  the  solution  to 


*  For  a  more  detailed  discussion  of  the  material  introduced  here,  cf.  ref. 
6,  p.  243,  and  ref.  7,  pp.  207-210. 
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Equation  (lid).  In  cartesian  coordinates,  the  substitution  of  (lie) 
into  (lid)  yields  the  result 


(  V  s)2]  +  [v  2  A  +  V(A  •  ^  ) 

A  V2S  +  2  (  V  A)  ‘  (V  S)  +  A(V  S) 


0 


(Ilf) 


The  latter  equation  will  be  satisfied  approximately  in  a  non-trivial 
fashion  if,  approximately,  the  following  equations  are  simultaneously 
satisfied. 


(VS)2  -  k 

V  2A  +  y  (A  •  ^  )  =»  0 

?  'w'U 

A  VS  +  2(VA)  •  (VS)  +  A(  \7  S)  *  k  ”  °* 


(gx) 

(g2)  (Hg) 
(g3) 


The  solution  S(x,y,z)  to  Equation  (Ilg^)  is  called  the  eikonal, 
and  Equation  (Ilg^)  is  called  the  eikonal  equation.  The  function  S  is 
important  in  optics  because  the  surfaces  S  =  constant  are  the  surfaces 
of  constant  phase,  and  they  may  be  thought  of  as  the  wave-fronts  determining 
the  mode  of  propagation  of  the  waves.  The  normals  to  the  surfaces  of 
constant  phase  determine  the  directions  of  the  "rays"  or  paths  along  which 
the  waves  propagate. 

Obviously,  if  k  and  A  are  constant  in  space,  the  system  (Ilg)  is 
satisfied  exactly  if 

(  VS)2  =  k  (hx) 

,  (Ilh) 

=  0.  (h2) 
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Now,  the  substitution  of  (He)  into  (lib)  yields  the  result 
Vk  , _ , 

A‘  k  +  ^7*  A  +  A  •  ^ S  -  0.  (Hi) 

Now,  if  for  anjr  reason  (e.g.,  a  constant  and  orthogonal  to  Vk)  we  have 
A  *  “k~  +  S7 '  A  "  °>  then  A  •  V  s  -  0,  which  says  that  the  field  direction 
is  transverse  to  the  direction  of  propagation,  since  VS  is  in  the  direction 
of  a  ray,  while  A  is  in  the  direction  of  the  field.  (The  case  in  which  k 
is  a  function  of  x  alone  and  A  =  ?yy  -  7  z  is  an  interesting  non-trivial 
example. ) 

Other  relations  analogous  to  the  above  can  be  derived  by  combining 
the  various  equations  in  other  ways.  The  assumptions  under  which  they  can 
be  derived  appear  to  be  far  less  stringent  than  the  assumptions  used,  for 
instance,  by  Sommerfeld  in  his  discussion  of  the  eikonal. 
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APPENDIX  III 


Reduction  of  Equation  (8) 

In  this  section  Equation  (8)  will  first  be  expressed  in  component 
form.  The  resulting  system  of  equations  will  then  be  transformed  to  another 
system  after  certain  restrictions  are  imposed  on  the  system. 

Let  u^,  u2,  u^  represent  the  three  coordinates  of  an  orthogonal 
coordinate  system  in  three-dimensional  space.  The  distance  ("differential 
arc  length")  ds  between  two  points  (u^,  u2>  u3)  and  (u^  +  du^,  u2  +  du2> 
u^  +  du^)  very  close  together  is  given  by  the  formula 
?  2  2  2  2  2  2 

ds  -  h^  du^  +  h2  du2  +  h^  du^  (Ilia) 

where  h^,  h^  h^  are  functions  of  u^,  u2,  u^  and  are  sometimes  called 
the  "scale  factors"  for  the  given  coordinate  system.  The  following  table 
1  ists  the  scale  factors  for  some  important  coordinate  systems  (cf.  ref. 


6,  p.  47  ff.): 

Coordinate  system  (u^,  u2,  u^) 

CM 

•C 

h3 

Cartesian  (x,y,z) 

1  1 

1 

Circular  cylindrical  (r,  0,  z) 

1  r 

1 

Spherical  (r,  4*  t  ®) 

1  r  sin  0 

r 

Elliptical  cylinder  (z,  ^  ^  ) 

IE 

!  c 

vl  -t{2  > 

1 12 

The  unconventional  ordering  of  the  latter 

two  coordinate  systems 

has  i 

used  here  for  later  reference.  In  order  to  express  the  components  of  (8) 
in  convenient  form  we  introduce  the  differential  operators  D2,  D3 
defined  by 
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} 


i  -  1,  2,  3, 


(IHb) 


where  <f>  is  an  arbitrary  function  of  u^,  u2,  u^.  Let  be  the  component 
of  E  in  the  direction  of  the  coordinate  u^  for  i  ■  1,  2,  3.  Then  (cf. 
ref.  6,  p.  50)  the  vector  equation  (8)  can  be  expressed  in  component  form 
as  follows 


From  Equation  (8)  it  also  follows  that  div(kE)  =  0,  which  can  be  written 

D1(h2h;jkE1)  +  D2(h3h1kE2)  +  D^h^kE^  =  0.  (IHd) 


In  the  light  of  the  remarks  made  in  Section  IIIA,  no  attempt  will 
be  made  here  to  reduce  the  system  (IIIc)  given  above;  but  we  shall  consider 
a  special  case  of  frequenc  occurrence  and  of  great  importantance, -namely;  the 
case  in  which  h^  =  1,  k  a  function  of  only,  and  the  ratio  of  h2  to  h3 
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is  a  product  of  functions,  each  of  which  is  a  function  of  only  one  of  the 
coordinates.  The  latter  assumption  is  satisfied  in  any  generalized 
cylindrical  coordinate  system  and  in  spherical  coordinates.  The  following 
useful  relation  can  then  be  established: 


Di(qDj)  <h3^>  -  iqDj  (qDi) 


for  i,  j  =  1,2,3,  i  *  j 
(IIIe1) 

where  (f>  is  an  arbitrary  function  of  the  coordinates.  The  proof  of  the 
latter  relation  proceeds  as  follows: 

<h3  i>  >  -  Di  ( s;  0  qh2^>-Di  ^  o-j 


’  Vi(M  )  +  hS  ir>  Wit >  +  o. 


'  iT  M  h?  D‘  )  0>2  t >  +  tt2  <#> ) 


q<DJ  ^ 

Di“j  (*'  S? 


(h2  (f>  ) 


By  virtue  of  the  assumption  about  the  ratio  h3/h2,  the  bracketed  quantity 
in  the  second  term  in  the  last  line  must  vanish,  yielding  the  desired 
result.  One  may  interchange  the  subscripts  2  and  3  in  (Ille^)  to  obtain 
another  relation,  namely, 


(h3  (j>  )  for  i,  j 


1*2,3 ,  i  4  j  • 

(IIIe2) 


By  suitable  manipulations  using  (Ille^)  and  (IIIe2),  one  may  prove  the 
following: 


(h3E3)  =  2(d1in  ^  )  D3(h2E3)  (fL) 


(IHf) 
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iz  n 

h3  1 


Ef  <h3E2>  -  °2  Di)  <h2E2>  ■  2<D1  q  <h3E2> 


Now,  starting  with  Equation  (IIIc^),  we  obtain 

h„ 


(f2) 

(Illf) 


°2  (  h23  D2j  +  °3  (  °3  )  +  h3h2k 


Bi  ’  °3  (  °1  )  <h3E3>  - 


Dj  ~ 


<h2E2>  “  r2  °1  (  Sf  °2  )  <h3E2>‘ 


Now  Equation  (Hid)  becomes,  with  h^  *  1,  and  k  a  function  of  u^  alone, 


t  OjO.jbjklj)  +  D2(h3E2)  +  D3(h2E3)  -  0. 

Solving  (Illg)  for  D^h^E,^)  and  substituting  the  result  into  the  last 
expression  in  the  last  line  above,  we  get 

°2  1  q'C)  +D3^D3l  *>3h2k 


(IHg) 


hT  °1  j  (W 


-  i^Dl(hJ  kDl)  (h2h3kEi )  V  hJDl)  1^D3  (h2E3)  ’ 


The  latter  equation,  when  rearranged,  yields,  with  the  help  of 


h2  °1  (  ^  £  Dx  J  (h2h3kE1)  + 


+  °3  f  r3  °3  )  ^3h2k 


q  hj  <h3E3>  •  q  Di  (q  -  2<Di  q  >  -MW- 


(IXIh^ 
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Similarly,  starting  in  a  different  way  from  (Illc^we  may  prove 


(h2h3kE1)  + 


h3h2k 


-  2  (Dl  J*n 


h3  *  °2  (h3E2)# 


(IHh2) 


Now,  in  cartesian  and  spherical  coordinates.  Equation  (Illh..)  or  (IIIh_) 

Ph2 

already  yields  an  equation  for  E,,  since,  in  these  cases  D,  (  X  n  —  )  = 

1  i  h„ 


0. 


In  circular  cylindrical  coordinates,  however,  a  complication  still  remains, 
h3  ^ 

since  -  D,  ^ n  —  *  D  ,/nr-^O.  In  order  to  avoid  even  more 
i-  h2  r  r 

tortuous  manipulations,  we  shall  at  this  point  make  the  additional  assumption 
that  all  variables  are  independent  of  This  is  tantamount  to  setting 

D2  =  0  wherever  T>2  appears  in  the  equations.  Equation  (IIIh2)  then  reduces 
to 


3  1 


h  DJ  I  (h2h3kE1)  +  (D3  D3  1  +  h3h2k)E1 


0. 


(IHj) 


Also,  we  obtain  immediately  from  Equation  (IIIc2>  the  equation 

,  \  h. 


+  D. 


h2h3  3 


(h2E2)  =  0, 


(Illk) 


an  equation  for  E2  alone.  Equations  (IIIc^)  and  (IIIc3)  reduce  to 
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By  further  algebraic  manipulations  the  latter  system  can  be  solved  to 


obtain  an  equation  for  alone.  From  a  practical  standpoint,  however, 
after  E^  is  obtained  from  Equation  (IIIj),  is  easily  obtainable  from 
Equation  (III>/^)  by  some  simple  integrations. 

We  present  here  the  elimination  of  E^  from  the  system  (111^)  for 
cartesian  and  circular  cylindrical  coordinates.  Then  h^  =  1,  and  h2  is 
independent  of  u^.  The  system  (III/  )  becomes 

(D32  +  k)Et  -  D3DlE3  (mL) 

(Him) 

Djti^Ej  =  (D1h2D1  +  h2k)E3.  (m2> 


Multiplying  (IHm^)  on  the  left  by  the  operator  D^h2D3,  (IIIm2)  on  the 

2 

left  by  the  operator  (D3  +  k),  subtracting  the  latter  from  the  former,  and 


simplifying,  we  obtain 


_  Dlh2Dl 


+  h„D_  +  h„k  E0  + 


(DLk) 


~  h2°3El  "  °* 


(IHn) 


Multiplying  (IHn)  on  the  left  by  the  operator  (D3  +  k),  and  using 

Equation  (Illm^),  we  obtain  finally 


(D32  +  k)  _D1h2D1  +  h2D32  +  h2kj  E 


(Dxk) 


3  +  — ~  h2D3  DLE3  =  0, 


(IIIp) 


thus  obtaining  an  equation  for  E3  alone. 
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APPENDIX  IV 


Special  Solutions  to  the  Wave  Equation 


In  this  special  case  in  which  the  electromagnetic  wave  is 
propagating  in  the  direction  of  the  gradient  of  k  (x-direction)  the  electric 
field  E  will  be  a  function  of  x  alone  and  will  satisfy  the  ordinary 
differential  equation 


+  kE  -  0. 


(IVa) 


When  k  has  the  function  form  given  by 


Ax  +  B 
(Cx  +  D)' 


;  A  -  AD  -  BC, 


(IVb) 


where  A,  B,  C,  and  D  are  constants,  then  Equation  (IVa)  has  the  solution 


C  1  Ax  +  B 

|DlAi  L’  ^2/3  Cx  +  D, 


+  D2Bi 


1  Ax 
^2/3  Cx 


«]] 


(CX  +  D), 


(IVc) 


where  D^  and  are  arbitrary  constants  and  Ai£^  Jj  and  Bi[^Jare  the  Airy 
functions  of  the  argument  ^  .  Now,  in  general,  it  is  not  possible  to 
identify  one  of  the  solutions  in  (IVc)  as  a  "forward-propagating"  wave, 
since  such  an  identification  would  require  an  examination  of  the  behavior 
of  the  functions  for  large  arguments.  However,  in  this  case  such  behavior 
need  not  be  relevant  to  the  problem  at  hand,  since  the  argument  may  not 
ever  become  large  in  the  region  of  interest  of  this  problem.  In  the  event 
that  k(x)  has  the  form  given  by  (IVb),  we  may  avoid  the  issue  by  discussing 


cf.  ref.  2,  pp.  44  and  45. 
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only  inhomogeneous  plane  sheets  of  finite  thickness,  bounded  by  semi¬ 
infinite,  homogeneous  media. 

The  following  example  is  discussed  in  Reference  3,  p.  229.  The 
example  is  originally  due  to  Schelkunoff. 

Suppose  E  to  be  given  as  a  function  of  x  of  the  form 

E(x)  =  cos  f3x  +  ee*^*,  (IVd) 

where  p  and  e  are  constants.  When  e  «  1,  this  function  describes  an 
essentially  standing  wave,  since  the  first  term  will  be  dominant.  However, 
this  function  can  be  written  in  the  form 

E(x)  =  ACxJe1  £  (X)  (e^ 

A(x)  =  1  +  e)2  cos2  px  +  e2  sin2  px  '  ( e 2)  (IVe) 

(j?  (x)  =  arc  tan  ^V'o  coV  px  *  (e3> 

The  latter  form  describes  a  "propagating"  wave. 
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APPENDIX  V 


Derivation  of  the  Reflection  and  Transmission  Coefficients 
_ for  an  Infinite.  Stratified  Slab _ 

Figure  1  shows  three  media.  Medium  I  extends  from  x  =  -  oo  to 

x  =  0  and  is  assumed  to  be  homogeneous  with  a  real  propagation  constant  k^. 

Medium  III  extends  from  x  =  d  to  x  =  co  and  has  a  propagation  constant. 

^3  =  ^3  +  *a3*  Medium  II  extends  from  x  =  0  to  x  =  d,  and  the  behavior 

of  the  electric  field  in  this  medium  is  governed  by  the  system  (9).  (As 

a  matter  of  fact,  the  behavior  of  the  electric  field  in  Media  I  and  III  is 

2  2 

also  governed  by  (9)  with  k  =  k1  and  k  =  k^  ,  respectively.)  The  X-Z 
plane  is  termed  the  plane  of  incidence. 

The  derivation  of  the  reflection  and  transmission  coefficients 
will  be  preceded  by  a  derivation  of  the  representations  for  the  reflected 
and  transmitted  waves.  Following  the  traditional  approach  (cf.,  ref.  6, 
pp.  492-494)  we  shall  consider  two  cases:  electric  field  vector  normal  to 
the  plane  of  incidence,  and  electric  field  vector  in  the  plane  of  incidence. 
The  reflection  and  transmission  coefficients  themselves  will  be  derived  from 
the  amplitudes  of  the  reflected  and  transmitted  waves  on  the  basis  of  the 
energy  propagation  relations  for  plane  waves  given  in  Reference  6,  p.  281. 

CASE  I.  Electric  Field  Polarized  Normal  to  the  of  Incidence 

_ (TE-Wave). _ 

In  this  case  E  =  E  =  0  in  all  three  media.  The  equation 

Z  X 

governing  the  behavior  of  E^  in  Medium  II  is  Equation  (9a).  As  indicated 

in  the  figure,  the  electric  field  in  Medium  I  consists  of  the  incident 

Ik 

wave,  having  the  functional  form  A^e  l(x  cos  0^  +  z  sin  0^),  and  the 
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ik 

reflected  wave,  having  the  functional  form  B^e  l(-x  cos  0^  +  z  sin  0^). 

The  electric  field  in  Medium  III  consists  of  an  outgoing  wave  (the  trans- 

ik 

mitted  wave)  having  the  functional  form  A^e  3  (x  cos  0^  +  z  sin  0^).  No 
attempt  will  be  made  at  this  point  actually  to  resolve  the  field  in  Medium 
II  into  forward  and  backward  propagating  waves.  Equation  (9a)  can  be 
solved  by  separating  variables.  If  the  solution  E^  is  represented  by 
Ey  =  ^  (x)e^YZ,  where  y  is  a  constant,  and  4>  (x)  is  independent  of  z,  then 


<j>  (x)  satisfies  the  ordinary  differential  equation 

+  (k  -  y2)  =  0. 

dx 


(Va) 


Let  4?  +(x)  +  B2  ^  _(x)  represent  the  general  solution  to  (Va), 

where  <f> +(x)  and  ^>_(x)  are  independent  solutions  to  (Va),  the  subscripts 
+  and  -  being  used  here  for  possible  future  identification  as  forward  and 
backward  propagating  waves.  Then  the  solution  to  (9a)  in  Medium  II 
can  be  expressed  as 

Ey  =  [A,  <f>+(x)  +  B2  <Mx)]  eiYZ.  (Vb) 

It  follows  from  Equations  (lb)  and  (3)  that 


S7  x  E  =  icu^H. 


(Vc) 


The  latter  equation,  when  expressed  in  component  form  in  the  case  we  are 


considering,  becomes 
dE 

-  =  i|j  o)H 

dz  'o  x 


(a) 


0  =  H 

y 

dE 

=  i^i  coH 
dx  o  z 


(b)  (Vd) 

(c) 
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Hence,  H  is  the  only  component  of  the  tangential  component  of  the 
z 

magnetic  field,  and  the  continuity  of  H  as  a  function  of  x  implies  the 
dE  Z 

continuity  of  as  a  function  of  x.  In  summary,  the  boundary  condi¬ 

tions  at  x  =  0  and  at  x  =  d  are 


E  is  a  continuous  function  of  x  at  x  =  0. 

y 


(eL) 


dE 

—1 

dx 


is  a  continuous  function  of  x  at  x  =  0. 


(e2> 


E  is  a  continuous  function  of  x  at  x  =  d. 

y 


(e3)  (Ve) 


dE 

_Z 

dx 


is  a  continuous  function  of  x  at  x  =  d. 


<v 


The  ultimate  goal  in  our  procedure  is  to  obtain  B^/A^  and  A^/A^,  from  which 
we  may  derive  the  reflection  and  transmission  coefficients. 

Application  of  conditions  (Vej)  and  (Ve2)  leads  to  the  following 
pair  of  equations: 

+  B1)eiklZ  8in  91  =  A2£+(0)  +  B24>_(0)J  eiYZ  (f^ 

(Vf) 

ikt  cos  01(A1  -  B^e^l28111  ®1  -[a^'cO)  +  B,,  <)>_ '  (0)J  eiYZ  (f2) 

where  the  primes  on  ^>+(x)  and<^_(x)  denote  differentiation  with  respect 
to  x.  Application  of  conditions  (V-e^)  and  (V-e^)  leads  to  the  pair: 


[^^(d)  +  B2f_(d)J 

[A2  ^j(d)  +  B2<f)jd)] 


iyz  _  .  ik_(z  sinQ,  +  d  cos  0q) 


(gL) 


(Vg) 


iyz  _ 


=  ik^cos  Q^A^e  3 


ik.(z  sin  0_  +  d  cos  0~) 


(g2) 
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When  we  multiply  through  Equation  (V-f^)  by  e  *^z ,  we  obtain  the  relation 
gi(kj.  8*n  ®1  "^z  =  constant.  The  latter  relation  must  hold  for  all 

values  of  z,  since  (V-f^)  is  to  hold  for  all  values  of  z.  This  situation 
can  prevail  only  if  sin  0^  -  y  =  0.  Similar  reasoning  on  (V-g^)  leads 
to  the  conclusion  that  k^  sin  0^  =  y.  In  summary, 

k^  sin  0^  =  y  =  k^  sin  0^.  (Vh) 

The  latter  can  be  regarded  as  a  generalization  of  Snell's  Law.  Also,  in 
view  of  the  relation  (V-h),  the  exponential  factors  can  be  cancelled  from 
(V-f)  and  (V-g),  yielding,  after  some  rearrangement,  the  following  set 
of  equations: 

a2  4>+(0)  +  b24>_(0)  -b1  -  ax 

I  t 

A2  4>+^  +  B2  +  ^1  cos  ®1B1  = 


(Vi) 

ik.dcos0 

A2^>+(d)  +  B2<£_(d)  -A3e  -  0 

t  ,  ik  dcos0- 

A2  <£+(d)  +  B2  -ik3cos  ®3A3e  =  0 

The  latter  set  of  equations  can  be  regarded  as  a  set  of  simultaneous  linear 

ik3d  cos  03 

algebraic  equations  in  the  unknowns  A2,  B 2,  B^,  A3e  .  The 

determinant  of  the  unknowns  is 


/N 

O 

V-/ 

+ 

+1 

0 

/-N 

O 

v-/ 

1 

-ik^cos0^ 

0 

i>+»> 

0 

+1 

(Vj) 

$+'<■» 

1) 

0 

+ik3cos03 
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*-x 

O 

X-' 

+ 

<t>_ '  (0) 

4+'<°) 

4>.'(o) 

i<0) 

<t>_(  0) 

♦+'<«> 

4>_'<d) 

+1^00802 

4>+<d> 

4>Jd) 

-ik^cos©^ 

<t>. '  (<*) 

-kjk^cos  S^cos  6^ 


<p+(  0)  <f>J  0) 

^+(d)  <p_(d) 


Let  the  minor  of  the  element  in  the  first  row  and  third  column 
of  A  be  denoted  by  D^,  and  let  the  minor  of  the  elements  in  the  second 
row  and  third  column  of  Abe  denoted  by  D£*  Then 


o 

/*N 

o 

X/ 

O 

X/ 

<P+'(  0)  t'(0) 

+^. 

/■N 

O 

1 

✓—v 

o 

x/ 

<fi+(d)  <{)_(d)  1 

“ik^cosQ^ 

- 

4>+  (d)  4>_  (d)  ik^cos  ©3 

<P+(d)  J>_(d)  . 

4>+’(d)  ^'(d) 

0^) 


(Vk) 


o 

o 

1 

o 

x^ 

<P+(0)  4>_(  0) 

4>+(  0)  4_(°> 

°2  = 

4>+(d)  (f>_  (d)  1 

=ik3cos@3 

- 

<P+  (d)  P.  (d)  ik3cos03 

</>+(d)  4>_(d) 

4+(d)  ^_’(d) 

(k2) 


Then,  (Vd)  A  =  +  (ik^cosB^jD^. 

Now,  solving  (Vd),  by  the  use  of  Cramer's  Rule,  we  find  that 
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-D^  +  (ikcos0^)D2 


B,  = 


1  +  (ikcosO^^  1 

2ik^cos0^  [>:  (d)  -  ik3  cos  03  4>_(d)] 

S  A1 

2ik^cos0^  [V  (d)  -  ik3cos©3  <£+(d)f 
- ^ -  Ai 


A2 


B2  = 


v(®1) 

V(m2) 

V(m3) 


<P+(0) 

4>_<o) 

-i 

1 

ik3dcos03  = 

■©- 

+ 

O 

N-/ 

4>_'(o) 

ik^cos0^ 

ik^cosO^ 

<fr+<d) 

0 

0 

4>+'  <d) 

4>_'(d) 

0 

0 

r- 

A 

2ik^cos0^ 

i - 

x-s 

CL 

v-/ 

£(d)  - 

<f>+'(  d)  4>_(d>J 

-ik  dcos0 

A  a  ^ 

A 

Ale 

V(m4) 

The  expression  |^>+(x)  <j>_  (x)  -  <£+  (x)  c />_(x)  J  is  the  Wronskian  of 

(j?+(x)  and  <p_(x)  an<*  therefore,  a  constant,  by  Abel's  Formula. 

CASE  II.  Electric  Field  in  the  Plane  of  Incidence.  (TM-Wave). 

In  this  case,  the  magnetic  field  is  normal  to  the  plane  of 
incidence,  =  0,  and  is  a  continuous  function  of  x  everywhere. 

Equation  1(a),  combined  with  (5)  and  (7),  yields 


V  x  H  =  — — - kE 

o*io 


*  cf.  Hildebrand,  Advanced  Calculus  for  Engineers,  Prentice-Hall, 
1949,  p.  31. 
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which,  when  expressed  in  component  form,  yields  the  following: 


aHy 

az 

SH 

3x 


i 


-i 


kE 


kE 


(r^) 


V(n) 


(n2) 


Since  E  is  a  continuous  function  of  x  for  all  values  of  x  and  z,  it 
z  1  a  H 

follows  from  V(n2)  that 


^  is  likewise.  In  summary,  we  may  say  that 


k  3 

X 

the  boundary  conditions  to  be  imposed  on  H  are  the  following: 

13  H  7 

^  are  continuous  functions  of  x  at  x  =  0  and  at  x  =  d.  In 


H  and 

y 


k  3 


addition,  we  note  that,  since  E^  is  a  solution  to  the  equation  9(b),  an 


equation  whose  variables  are  separable,  we  may  express  E  in  the  form 


Ex  =  /<x)  e1  % 


V(0) 


where  ^  is  a  constant,  andy~(x)  is  a  function  of  x  alone  satisfying  the 
differential  equation 


-k  +  o. 


V(p) 


Combining  equations  V(n^)  and  V(0),  we  get 


3H  i  t  : 

jf  =  Z^k^x)  e  * 


which,  upon  integration  with  respect  to  z,  yields 

H  =  — - - -  k^(x)  e1  ^  Z  =/?(x) 

y  £  a*i0  r  '  ' 


where  we  have  introduced/^  (x)  for  k  ^(x) 
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The  procedure  in  this  case  then  exactly  parallels  the  procedure 


« 


1 


used  in  Case  I,  except  that  we  are  now  talking  about  magnetic  fields 

instead  of  electric  fields.  In  place  of  the  symbols  Ay  A y  Ay  By  By 

<£+(x),  <pjx),  A,  Dy  Dy  use  the  symbols  Oy  Oy  ay  0y  P2, /5+(s), 

(x),  §  ,  Gy  Gy  respectively,  where  the  a's  and  p's  are  now  magnetic 

field  amplitudes.  Because  of  the  slight  difference  in  the  form  of  the 

boundary  conditions,  it  turns  out  that  a  further  modification  arises  in 

that,  in  all  of  the  determinants  in  question  (but  not  in  exponential 

expressions)  should  be  replaced  by  — and  k^  should  be  replaced 

by  Thus,  we  obtain,  for  example,  comparing  V(m,  )  and  V(j), 

k3  4 


2i 


M2i 


cos  0, 


a3  =  - 


/%(d)  /f.  (d)  -/*+  (d)  (d)| 

: - 5 - 


a^e 


■ik^dcosS^ 


V(r) 


4<o) 

/t<0) 

1 

0 

4  (°) 

A '  <°> 

-ik(O)cos0^ 

A 

kl 

u 

4(d) 

/?(  d) 

0 

1 

4<d) 

' (d) 

n 

ik(d)cos@2 

u 

k3 

V(s) 


Reflection  and  transmission  coefficients.  Following  the  methods  used  in 
Ref.  6,  p.  496,  we  have  the  reflection  and  transmission  coefficients  R1 
Ty  ,  Ry,  Ty  for  cases  I  and  II  given,  respectively,  by  V(t)  and  V(u) 
as  follows  (assuming  k^  real): 
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'X 


ci 


-Dj.  +  (ik1cos81)D2 
Dj  +  (ik1cos01)D2 


ReCk^cosS^) 
k^  cos  0^ 


(tl> 


V(t) 


(t2) 


4k1cos01Re(k3cos0^) 

I  4k  d  cos  9  |  2  f*W  *'W  ■  V <d)  *-(d)  I 


11 


L11 


£i 

2 

"Gj.  +  (ik(0)/ki  cos  0^)G2 

a 

P 

Gj^  +  (ikCOJ/kj^  cos  91)G2 

(u,  ) 


V(u) 


a. 


a. 


CO80, 
COS  0. 


O  CM 
^  i-M 

/Vd)  -  f+W  /^.(d) 

2 

S  gk^d  COS03 

1  2 

0080^0080^ 

(U2> 
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